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I. INTRODUCTION 

A central problem in making Standard Model predictions for K meson decays, including 
the important CP violating amplitudes, is calculation of the decay K — > tttt into the 1 = 
final state of two pions. In a standard calculation in which the K meson is at rest, the 
contribution of the two pion final state is expected to be very difficult to extract from a 
lattice QCD calculation. The physical state has the same quantum numbers as the state 
with the two pions at rest (the difficulty emphasized by Maiani and Testa Of even 
greater concern is the fact that the QCD vacuum also has these same quantum numbers. A 
very appealing approach to deal with these problems is to study the decay of a K meson 
with non-zero momentum. For a calculation in finite volume, it is possible to adjust the 
momentum of the K meson and the box size so that both the transition amplitude to the 
vacuum must vanish (since the vacuum has zero momentum) and the final two-pion state 
must have physical relative momentum and total energy equal to that of the decaying K 
meson. 

For such a calculation to be useful we must be able to relate the decay amplitudes com- 
puted in finite volume using this technique to the infinite volume matrix elements determined 
by experiment. To a large extent the effects of finite volume on the particles involved are 
relatively mild. Typically the volumes used in such a calculation can be chosen sufficiently 
large that they do not significantly distort the K and 7r particles whose physical size will be 
much smaller than that of the box employed. The most obvious finite-volume effect will be 
the quantization of the energy levels of the two-pion final state. This is actually an advan- 
tage, permitting the box size to be adjusted to make one of the discrete, two-pion energies 
match precisely that of the K meson. For the case at hand, this can even be the lowest 
energy tt — tt state. 

However, the violation of rotational symmetry by the finite-volume boundary conditions 
does induce an important distortion in the computed decay rate. Because of the resulting 
non-conservation of angular momentum, the two-particle state into which the decay occurs 
is actually a mixture of states with many angular momenta. For typical lattice volumes 
the actual decay of the K meson into these higher angular momentum states will be very 
small — angular momentum is effectively conserved at the short distances over which the 
decay occurs. However, the presence of these extra angular momentum states effects the 
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normalization of the physical J = amplitude which appears in the matrix element. 

This finite volume normalization problem has been solved by Lellouch and Luscher [2] for 
the case of the decay of a K meson at rest. This problem has also been solved by a different 
approach in Ref . p] • In this paper we will generalize the approach of Lellouch and Luscher 
to obtain a result for states with non-zero total momentum. Central to their argument is an 
earli. t— of Lilscher flfl which defines the allowed, finite-volume, tw„-pa rt ic.e 
energy eigenvalues in terms of the infinite volume, two-particle scattering phase shifts for 
energies below all inelastic thresholds. This discussion must also be generalized to the case 
of non-zero center-of-mass momentum. 

This topic has been studied earlier by Rummukainen and Gottlieb jfj. Their treatment 
involves an application of relativistic two-particle quantum mechanics. We believe that it is 
also important to study this problem starting from the equations of quantum field theory. 

Our strategy follows closely that of Luscher 0, 0] and Lellouch and Luscher We 
first discuss the energy quantization of finite-volume, interacting, two-particle states with 
non-zero center-of-mass momentum. Following Luscher, this is first done in standard, two- 
particle, non-relativistic quantum mechanics in Sec.|n] In Sec. IHII we begin with the Bethe- 
Salpeter equation of relativistic field theory and, again following Luscher, show how this 
equation when restricted to a particular 7-dimensional subspace of the 8- dimensional, 2- 
particle momentum space reduces to the standard Lippmann-Schwinger equation describing 
the earlier non-relativistic system. Finally, in Sec IIVI we use this result to generalize the 
argument of Lellouch and Luscher to the determine the finite volume corrections to the 
decay amplitude computed for states with non-zero total momentum. 

The issues addressed in this paper have also been considered by Kim, Sachrajda and 
Sharpe. Using a related but different approach, they have also confirmed the validity of the 
results of Ref. ly] and derived the generalization of the result of Lellouch and Luscher for 

Q 

the case of non-zero total momentum. Their paper p| is being released simultaneously with 
the present article. 

II. FINITE VOLUME, NON-RELATIVISTIC, TWO-PARTICLE STATES 



We begin by considering a simple, non-relativistic system of two distinguishable particles 
confined in a cubic box of side L and obeying periodic boundary conditions. The system 
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is described by a wave function ip(ri,r 2 ) which is periodic in r x and r 2 separately. An 
eigenstate of energy ip E obeys the Schrodinger equation: 

{~ - |f + V(\n - f 2 |)}^(f 1? f 2 ) = EMn,r 2 ), (l) 

where m is the identical mass of the two particles and V{\f\ — r*2 1 ) their rotationally invariant 
interaction potential. 

Such an equation is conventionally simplified by changing to center-of-mass and relative 
coordinates: 

R = (2) 
f = r x -r 2 (3) 

with conjugate momenta 

P = Pi + P2 (4) 

P1-P2 f . 
P = —3—, (5) 

where pi is the momentum conjugate to the coordinate fj. 

Using ipg e ^ (r, R) to represent the original wave function expressed in terms of f and R, 
we can write down the standard equation which it obeys: 

{- - S + y(|rlM^M) = (6) 
L 4m m J 

The periodicity of ^(fi,^)^'' under the simultaneous translation r. \ — > f! j + e^L (where 
is a unit vector parallel to one of the edges of the box) implies the periodicity of the 
wave function ip r under a translation of R by L. Thus, the conserved total momentum P 
must obey the quantization condition: P = J2 k=1 2nnkek/L for integer n^. It is easy to see 
that in a direction k for which the integer n& is even, the corresponding component of the 
relative momentum p k = 27in' k /L while if is odd then p k = 27r(n' fc + \)/ L, where n' k is an 
integer. With this change of coordinates and a specific choice of P, our two-particle problem 
reduces to the quantum mechanics of a single particle in an L 3 box obeying either periodic 
or antiperiodic boundary conditions on each of its three opposing faces. 

In infinite volume, this can be viewed as a scattering problem often phrased as a 
Lippmann-Schwinger integral equation. One defines an energy eigenstate tp 1 ? whose incom- 
ing part (that term with radial dependence e~ ipr ) is that of a plane wave with momentum 
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p. Necessarily, the outgoing part of ip 1 ^ is more general, being created by scattering from 
the potential V. Such a state must have energy E cm = p 2 /m. Manipulation of Eq. El easily 
produces the desired integral equation: 

where (j)p is the plane wave solution <fip{f) = e l ^' T of the free Schrodinger equation and 
Hq = — V 2 /m, the free Hamiltonian of a particle with the "reduced mass" m/2. Of course 
the full solution to Eq. ip^ el \f, R), is a product of a plane- wave depending on the center- 
of-mass coordinate R and the wave function above: 

4 e \r,R) = e ip - M ^(?) (8) 

and the total energy E is related to the energy in the center-of-mass system by E = E cm + 
P 2 /4m. 

Examining the asymptotic behavior of Eq. [7J one derives the standard relation between 
the conventional scattering amplitude f(9) and the matrix element of V between the plane 
wave state 4>p> an d ip l p'- 

f(6) = -2n 2 m(<J ) p,\V\^) J (9) 

where p' • p — p 2 cos(6 ) ). We obtain an equation closer to the relativistic Bethe-Salpeter 
equation by defining the T matrix as: 

(p'\m = (<j ) p,\v\ijf). (io) 

Here we are using the conventional Dirac bra-ket notation to represent momentum eigen- 
states: (f>p= \p). Note, Eq. [TUl defines matrix elements of T even when \p'\ ^ \p\ and energy 
is not conserved. 

Using the matrix T we can rewrite Eq. m a form that will be useful later, if we multiply 
by V and transform to momentum space: 

/d^k -> 1 

This equation explicitly involves the matrix elements of T between states with different 
energies. Equation is a component of standard scattering theory and applies only to 
the case of infinite volume. Since we are interested also in the eigenfunctions and energies 
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for the finite volume problem it will be helpful to cast the standard Schrodinger equation, 
obeyed by an eigenstate, ip n (r) with the discrete energy 

{H + V)^ n = E™i> n , (12) 

into a similar form: 

n 

/d 3 k -> 1 

This equation demonstrates that the state VipE n solves the homogenous Lippmann- 
Schwinger equation if the energy argument in the denominator, E cm = p 2 jm is continued 
to that state's actual energy E cm = E™. Equation HH1 can be used in infinite volume to 
determine the energy of a possible bound state. It also can be applied to the case of finite 
volume if one makes a simple replacement of the integral over the relative momentum k by 
an appropriate discrete sum. 

In Eq. m we have followed the standard procedure, exploiting the separation of relative 
and center-of-mass variables permitted by Eq. |Hl and written that integral equation as 
an equation obeyed by functions of a single, three-dimensional, relative momentum k or 
p. The three-momentum of the center-of-mass, P, disappears from the problem once we 
use the energy in center-of-mass system, E cm . However, we could create a more explicit 
analogy with the relativistic discussion to follow if we viewed the states \p) and \k) as 
functions of the four-momentum of the center-of-mass as well: P(p) = (£- + p 2 / m, P) and 
— (j~ + k 2 /m, P) respectively. Had we done this, Eq. [TT] would be an equation obeyed 
by functions of seven variables. Each factor in this equation would be explicitly diagonal 
in P and the difference of these total energy variables Pq(p) and Po(k) would replace the 
present denominator in that equation. (A similar remark applies to Eq. EI] as well.) 

The final step to be reviewed in this section is the connection between the infinite volume 
scattering problem, defined by Eq. ^2 an d the finite volume energy eigenvalues of the original 
Schrodinger equation, Eq. El or Eq. UHl This is the problem solved by Luscher in Refs. 0,0]. 
Recall that the scattering amplitude f(6) can be written as a sum over partial waves as: 



fie) = J2( 21 + l ) e -^- p ii cos (P)i ( 14 ) 



1=0 * 
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where the Si are the standard scattering phase shifts and the Pi(cos(9)) the usual Legendre 



polynomials. In Refs. 



2Q, 



Liischer examines the case of a potential of finite range, V(r) = 
for r > -R max , and for the case L > 2_R max derives a relation between the allowed energies 
E cm = p 2 jm in the finite box and the phase shifts b\. For the simplest case where all b\ ~ 
for / > 0, this finite-volume quantization condition becomes: 

nir - 5 (p) = (j)(q) (15) 



where n is an integer, p = \JrnE, q = pL/2ir and the function <fi(q) is a known kinematic 
function given by: 

3 /2 

tan0(g) = - J 9 0(0) = 0, (16) 

^oo(i;r) 

with the zeta function Z 00 (s; q 2 ) defined by 

Z 00 (s;q 2 ) = ^=J2( n2 - ( i 2 y S - 

* / A Tr ' ™ 
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The zeta function defined above applies to the case that the integers appearing in the center 
of mass momentum P are even. For the case that one or more is odd, the finite volume 
problem will obey anti-periodic boundary conditions in those directions and an appropriate 
offset of 1/2 must be added in the summation in Eq. E3 as discussed in Ref. Jf|. We 
conclude that the total energy of a finite volume system with total momenta P, is given by 
E = P 2 /Am + E cm where p = VmE cm obeys Eq. EH1 

Thus, Luscher's relation for non-relativistic quantum mechanics between the scattering 
phase shifts and the energy eigenvalues for that same system in a finite volume is straight 
forward to generalize to the case that the two-particle system carries non-zero total momen- 
tum. We now consider the generalization of the next step to non-zero total momentum: 
the connection between the Bethe-Salpeter equation of relativistic quantum field theory and 
Eq. HU above, obeyed by the non-relativistic T matrix. This represents the new result of 
this paper. 



III. FINITE VOLUME RELATIVISTIC TWO-PARTICLE STATES 

In this section we will generalize to the case of non-zero total momentum the procedure 
introduced by Liischer to reduce the Bethe-Salpeter equation of relativistic field theory to an 
equation whose form is identical to the non-relativistic Eqs. and ED This will permit the 
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quantization condition described in Eqs. [ToT[T7| to be applied to quantum field theory, and 
QCD in particular. As discussed in Refs. [J, we must do this by carefully distinguishing 
effects of finite volume which fall as powers or exponentially in the system size. Note: our 
objective is to derive an equation which is both similar in form to Eqs. ^2 and and also 
accurate for both finite and infinite volume so that it can be used to relate the finite volume 
spectrum and the infinite volume scattering amplitude. 

We begin with the standard Bethe-Salpeter equation, which connects the amputated 
four-point function, T(p' 1 ,p' 2 ;pi,p 2 ), the two-particle irreducible kernel K^p'^p'^px^p^) and 
the single particle propagator, A(k 2 ): 

T(p[,p' 2 ;Pi,P2) = K(p' 1 ,p' 2 ;p l ,p 2 ) 

+ / -^ K (Pi,P2;P/2 + kP/2-~k) 

■A((P/2 + fc) 2 )A((P/2 - k) 2 )T(P/2 + k, P/2 - k, P i,p 2 ), (18) 

where P = pi+p 2 = Pi+P 2 is the total four-momentum. For simplicity, we have also removed 
an overall ^-function for the conservation of total four-momentum from T(p' 1 ,p 2 ;pi,P2), 
Kij/iiPh'iPiiPz) an d A(A; 2 ). For example, T(p' 1 ,p' 2 ;pi,p2) is defined by: 

(27r) 4 5 4 (p' 1 +p 2 -pi -P2)T(p[,p' 2 ;pi,p 2 ) = 




■(Ol^x'^ix'^ixMx^lO)^ (19) 

where 4>{x) is the quantized scalar field used in this example (normalized so that the single 
particle pole in the 2-point function has unit residue), (. . .) con n indicates that only connected 
diagrams are to be included and we are following the conventions of Peskin and Schoeder j^] . 

For simplicity we write the sum over the internal four-momentum k in Eq. as an 
integral. However, this equation applies equally well to a finite volume system if the spatial 
part of this continuous integral is replaced by an appropriate discrete sum. Specifically, 
in finite volume the total momentum operator P op is conserved and, given periodic spatial 
boundary conditions, has components which are quantized in units of 2tt/L. The single- 
particle three momenta, p, pi', P/2 + k, etc. correspond to indices in a finite-volume 
Fourier series and are also quantized in units of 2n/L. 

Equation^Jis more general than needed, constraining the off-shell, two-particle scattering 
amplitude, a matrix acting on a space of functions of eight momentum components. We 
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must specialize this equation, reducing the number of momentum variables from eight to 
seven. The choice of this seven-dimensional restriction of Eq. El is the key step in the 
desired generalization to non-zero total 3-momentum, P ^ 0. Once that has been done, 
the resulting equation will still be quite different from our non-relativistic target, Eq. ITT1 
However, the final steps connecting Eq. El with Eqs. ^2 and El proceed in a fashion very 
similar to those in Luscher's original derivation: We extract from the second term in Eq. El 
particular pieces which are regular functions of the total energy P . By a re-arrangement 
procedure, these terms are incorporated in a modified kernel K(p' 1 ,p' 2 ;pi,p2)- After this 
step, our re-arranged equation will have the same form as the non-relativistic Eqs. E] and 

o 

In contrast to the non-relativistic case, our new equation will have a volume dependent 
potential. However, if we remain below the four-particle threshold and have introduced into 
our re-arranged kernel K only regular factors, the difference between the finite and infinite 
volume kernels will vanish exponentially in the box size. Similar exponentially small errors 
will come from the failure of the resulting potential to have a truly finite range. Thus, we 
will be assured that the quantization condition in Eq. El will apply to the relativistic case 
with P 7^ up to exponentially small corrections. 

We begin by restricting the general Bethe-Salpeter equation of Eq. El to a carefully 
defined seven- dimensional momentum surface, making it similar to the target Eq. EJ (As 
discussed earlier, Eq. ^2 can also be viewed as involving functions of seven dimensions 
provided the trivial dependence on the total four- momentum is included.) 

First express the pairs of two particle momenta p[,p 2 and p%,p2 using relative and total 
momenta: 

pi )2 = P/2 ± k 

p' 1)2 = P/2 ± k'. (20) 

In Eq. El we have already imposed four- momentum conservation, removing an over-all mo- 
mentum conserving delta function from both T and K. Next we impose a further condition 
on the four- momenta k and k! requiring that ko = ftkn, where (3 = \P\/Pq and k\\ is the 
spatial component of the four-vector k in the direction of P. A similar condition defines 
a restricted value for k! and the integration variable k in Eq. El This condition is not 
immediately useful since, while we can consistently impose it on the external k and k' vari- 
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ables in Eq. the integral (or sum) over k does not obey any such restriction. Note, the 
restriction k = (3k\\ can be applied equally well in finite or infinite volume since, while k\\ 
will be discrete for the finite volume case, the time component, k , is always continuous and 
can be chosen to obey such a relation. 

To make progress we must remove from the integral over k some terms which are not 
singular as Pq approaches an allowed two-particle energy, P — > uo + + where the single 



particle energies are given by u± = y (P/2 ± k) 2 + m 2 . This can be done by generalizing a 
discussion in Ref. j^J and arguing that the portion of the product A((P/2 + A;) 2 )A((P/2 — k) 2 ) 
which is a singular function of Pq comes from the product of the single-particle singularities 
in each of the A factors. This singular term, arising when these poles pinch the k contour, 
can be viewed as a distribution in k , allowing us to write this product as: 

— 1 IT 

A((P/2 + *) 2 )A((P/2-*) 2 ) = P 2 /4 + k ,_ m 2 + l J ( p - k ) + mk) 

= S(P,k)+R(P,k), (21) 

where the function R(P, k) is a regular function of Pq in the interval below the four-pion 
threshold: 

fp2 jp2 

2W— +m 2 <P <4W— + m 2 (22) 

and the singular part, S(P, k), is defined by Eq. |^ Equation |^ is derived in the Appendix 
and is a generalization of Luscher's Eq. 3.16 from Ref. jj] to the case of P 7^ 0. 

The next step absorbs the contribution from the regular function R(P, k) as follows. First 
rewrite Eq. i n a more symbolic form exploiting the decomposition in Eq. 

T = K + K(S + R)T. (23) 

The second term on the right-hand-side can then be moved to the left hand side: 

(1 - KR)T = K + KST. (24) 

Finally we divide by the factor (1 — KR) and define 

K = K (25) 

l-KR y ! 

The resulting equation can then be written: 

T(k'- k) = K(k'; k)+f ^uK{k>, k) ~r b f - P) T{k, k). (26) 
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Here we have replaced the variables Pi and p\ with the total and relative four-momenta P, 
k! and k and suppressed the variable P. 

It is now easy to see that Eq. ESI has a form identical to the original non-relativistic 
Lippmann-Schwinger equation, Eq. ^2 First we observe that the delta function 5(k ■ P) 
forces the integration four-momentum k to obey our restriction: 

= k ■ P = k Po ~ k\\\P\ or k = pk ll . (27) 

This permits us to impose this relation between the time and parallel components of the rela- 
tive four-momenta everywhere in this equation, effectively reducing it to a three-dimensional 
integral equation as is the case for the non-relativistic problem. (As in that case, this equa- 
tion is diagonal in the total four momentum, P, the remaining four of our seven-dimensional 
momentum variables.) 

Second, we observe that the denominator has a non-relativistic form if we rescale the axis 



parallel to P by a factor of 7 = 1/ \A — (3 2 : 



p2 p2 p2 

— + k 2 -m 2 = -p + kl -% 2 - m 

4 4 4 u 



p2 p2 I 

T-T-^- 1 "-'" 2 - (28) 



Thus, if we change variables from fcii to 



h = ~ k h (29) 
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the denominator has the normal Laplacian form. We obtain a complete match between the 
denominators in Eqs. ^2 and EH and that in Eq. ESI if we remove a factor of m from the 
denominator of Eq. ESl an d identify the non-relativistic energy E cm with (P 2 — P 2 )/4m — m. 

Here we should recall the standard connection between Bethe-Salpeter equation, Eq. 
and the version of the Schrodinger equation given by the homogenous Eq. EH which deter- 
mines the discrete, finite- volume energies. Since the Bethe-Salpeter equation, e.g. Eq. ESI 
holds in finite volume, the energy eigenvalues, E n for the finite- volume, relativistic system 
will correspond to poles in the 2-particle scattering amplitude T obeying that equation. 
As one approaches the singularity at P ~~ > E n the inhomogeneous term, proportional to 
the kernel K, is not singular and can be dropped from the equation leaving a homogenous 
equation identical in form to Eq. E^ 
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Finally we must investigate the rotational symmetry of the kernel, K(k'; k). Since this 
function also depends on the four-vector P, there is possible rotationally asymmetric depen- 
dence on arguments of the form P ■ k' and P ■ k in addition to the acceptable dependence on 
(k') 2 , k 2 and k! ■ k. Fortunately, if the components of k! and k parallel to P are rescaled as 
described in Eq. |2HJ it is easy to see that the resulting function K(k'; k) becomes rotationally 
symmetric. This can be demonstrated by exploiting the Lorentz invariance of the function 
K(k'; k) which is assured by the covariant separation of regular and singular parts used in 
Eq. EU We use this Lorentz symmetry to equate K{k'; k) to its value in the rest system: 

K(P\ k ; k) = K(P cm ; k cm \ k cm ), (30) 

where for clarity we now also display the dependence on the total four-momentum P. 

Since in the center of mass system P cm = 0, the right hand side of Eq. Ellis a manifestly 
rotationally invariant function of the center of mass variables k' cm and k cm . However we can 
easily express these variable in terms of the original k' and k. For example: 

— * — * 

(k cm )_L = k± (31) 

(k cm )\\ = -Pko) = -k\\ (32) 

(k cm )o = 7(^o - pk\\) = 0, (33) 

where the second and third lines follow from the constraint obeyed by ko, imposed when we 
specialized to this three-dimensional equation. The third equation for (fc cm )o is necessary to 
insure that rotational asymmetry is not introduced through a non-symmetric dependence of 
this variable on the original laboratory variables. 

Note that the rescaled variables on which K(P; k'; k) depends symmetrically, (fccm)|| = 
(1/7)A;|| and {k' cm )\\ = (1/7)^11, are precisely the variables in terms of which the denominator 
in Eq. [2H1 contains the standard Laplacian. 

For clarity, we now rewrite the resulting three-dimensional integral equation, Eq. |2HJ in 
terms of these translated momentum variables, labeled suggestively k CUL : 

(^cml-^l^cm) (^cml-^'l^'cm} "H / ~7n ( ^ C m I ^ I ^cm ) 777^ 5 (^cm | T\ k cm ) . 

J (2ny 2P 0(p2/ 4 _ m 2_ A ,2 m ) 

(34) 

> — # . 

While the variables k' cm , k cm and k cm may appear to be Lorentz transformed, center-of-mass 
variables, in fact, they are simply the original variables defined in the laboratory system 



13 



except for a transformation of scale for the single component parallel to P. This rescaling is 
equally well defined if the corresponding variable is continuous or discrete. The only explicit 
use of Lorentz symmetry is to constrain the possible dependence of the function K(P, k', k) 
on its arguments. 

Thus, we have demonstrated that our original, field-theoretic Bethe-Salpeter equation 
can be rewritten as a non-relativistic Lippmann-Schwinger equation if we change to rescaled 
variables {k cm )\\ = (1/7)^11. We must now ask if such a rescaled set of momenta corresponds 
to an actual finite volume problem. Recall that we began by examining a two-particle 
problem in a finite, cubic box of side L. The transformation to the coordinates R and r of 
Eqs. El and El and the choice of total momenta P then requires that the relative momenta k 
or k' that appear in the Bethe-Salpeter Eq. |2""]have the form 2ir(nx, n 2, TI3) / L where rij is an 
integer or half-integer depending on whether LPi/2n is an even or odd integer. Now, we have 
recognized that this problem is equivalent to a non-relativistic problem with transformed 
momenta given by Eqs. H"""] and [221 If these k cm momenta correspond to those for a finite 
volume problem, then we have succeeded in casting the original relativistic problem into a 
non-relativistic problem which can be solved using the techniques developed by Luscher in 
Refs. QQ. 

That this is in fact the case can be seen by examining two simple examples. In the first 
example P = 2tc(0, 0, 1)/L. In this case the transformed variables fc™ 1 take on a very simple 
form: 

ki m = k x = (35) 

kT = k 2 = ^n 2 (36) 

kT = 1 (k 3 -pk ) = -k 3 = ^n 3 (37) 

7 7^ 

where ni, n 2 and + ~ are integers. These quantized momenta correspond to a simple finite 
volume of length L in the 1-and 2-directions and expanded length 7L in the 3-direction. If 
periodic boundary conditions are imposed in the 1- and 2-directions and anti-periodic condi- 
tions imposed in the 3-direction, the resulting quantized momenta will correspond precisely 
with those in Eqs. 133113*71 Thus, after generalizing Liischer's non-relativistic technique to 
this sort of asymmetric box, we will obtain the desired relation between the infinite-volume 
scattering phase shifts and the discrete finite volume energies in this asymmetric box. Of 
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course, our analysis above has determined that these same energies will be found in the 
original relativistic problem with total momentum P = 27r(0,0, 1)/L and cubic box of side 
L. This agrees with the result obtained Rummukainen and Gottlieb |6] and corresponds to 
of immediate practical interest. 
Now let us examine a second case where P = (1, 1,0)2tt/L. Using Eqs. RTTI and PT21 the 
rescaled momentum k cm is given by 

_ _ k ■ P _, 1 k ■ P _ 

k ^ = k-^-J-P + -^J-P ( 38) 

|p|2 7 |p|2 

where the first two terms correspond to the untransformed perpendicular component while 
the third term is the transformed parallel piece. Written in terms of the individual compo- 
nents this equation becomes: 

27T 

fc 3 cm = ^3 = ^s T - (41) 

These quantization conditions can be easily realized if we impose the following 
(anti-)periodicity conditions on the wave function ip cm (f) on which the operators in the 
integral equation, Eq. OH act: 

f m (r] = -V> cm (f+A), for i = l,2 (42) 
f m (rl = +4j cm {f+e 3 L) (43) 

where the displacement vector Dj is chosen to pick out the integer rii from the dot product 
Di ■ k cm = 2nniL: 

^i = e!^-L + e 2 ^-L (44) 
J D 2 = e 1 ^L + e 2 ^L. (45) 

The quantization condition given in Eqs. 021 and HH1 is equivalent to requiring that the 
wavefunction -0 cm (r) obey anti-periodic boundary conditions on the faces of a rhombus whose 
sides are parallel to the vectors Di and whose diagonals have length | D\ + Do\ = \/2jL and 
\D± — Z?2 1 = V%L, precisely the earlier result of Rummukainen and Gottlieb p]. From these 



15 



two examples, it is clear that the case of general total momentum P = (ni, n.2, n^)2T[ / L can 
also be realized by imposing (anti-)periodic on an appropriately distorted volume. 

The results of this section can be summarized by returning to our first example, the case of 
a symmetrical L 3 box and a two-pion state with total momentum oriented in the 3-direction: 
P = (2n/L)e3. Under these circumstances, the Bethe-Salpeter equation obeyed by the two- 
particle, off-shell scattering amplitude has been shown to be equivalent to a Schrodinger-like 
wave equation obeyed in an asymmetrical box with sides L x L x 7L with the longer 7L 
side parallel to the 3-direction. Thus, up to exponentially small corrections, the energy 
eigenvalues of the original 2-particle, L 3 system can be predicted using this L x L x jL, 
Schrodinger-like system. As a result, the allowed energies E of this original system must 
obey a quantization condition similar to that given in Eq. I15lwhere the function 4>(k) in that 
equation must be modified to describe the anti-periodic boundary conditions and expanded 
length in the 3-direction: 

nn - 5 (k cm ) = (f>(q) (46) 

f m( \ l^ l2 q 
tan 0(g) = 



Z 00 (l;g 2 ;7)' 

0(0) = (47) 
^oo(s; q 2 ; 7) = -4=E( n ? + ^+ ~ 2 (n 3 + \f - q 2 )-\ (48) 



where n is an integer, k cm = h^/P 2 — (2m) 2 and q = k cm L/2ir. This is the original result of 
Rummukainen and Gottlieb 



IV. TWO-PARTICLE DECAY OF STATES WITH NON-ZERO MOMENTUM 

The last part of this discussion is a generalization of the arguments of Lellouch and 
Liischer in Ref. |2j to the case of non-zero total momenta. Fortunately, this is very straight 
forward because the methods employed in that paper work equally well for P ^ once 
the formula relating energy levels and scattering phase shifts has been generalized to this 
non-zero momentum case. 

We begin by reviewing the Lellouch-Luscher approach, which is somewhat indirect. For 
our present purposes we will describe this method for the case of non-zero center-of-mass 
momentum, P ^ 0. One considers a finite volume system with both a K° meson and a 



16 



degenerate 7r + 7r state, each with total four-momentum P. Here the finite volume has been 



adjusted to insure that E nn = \J m 2 K + P 2 , including the effects of the tt — tt interaction. 
Next, the effects of the weak interaction Hamiltonian, Hw, mixing these two states, are then 
examined in perturbation theory. To zeroth order, the K° and tt + tt~ states are degenerate 
and non-interacting. To first order these states mix and their energies can be determined in 
first order, degenerate perturbation theory: 



P o = \Jm 2 K + P 2 (49) 

-> Po + AP (50) 

= P q ±(ti + 7i-\H w \K } (51) 

= P ±M, (52) 

where the states appearing in this formula are finite- volume states with non-zero total mo- 
mentum, normalized to unit probability. 

The next step relates the finite volume amplitude M, which can be computed directly in 
a lattice QCD calculation, with the infinite volume matrix element of Hw that determines 
the physical partial width. This step uses Eq. 14*61 to relate the infinite volume phase shift, 
computed at the quantized, finite-volume energy P + AP with that energy shift related to 
M by Eqs. 15*0115*21 Since the effect of H w on the scattering phase shift can be determined 
analytically in terms of the infinite- volume matrix elements of Hw, this equation will then 
relate the known, finite volume matrix element of Hw with the desired, infinite volume 
matrix element. 

Thus, we must compute the variation in the tt — tt scattering phase shift caused by the 

resonant scattering into the K meson state. This infinite volume calculation is most easily 

done in the tt — tt center of mass system and follows easily from the single tt — tt scattering 

diagram with a K meson intermediate state shown in Fig. [*Q giving: 

k cm \A\ 2 
327rm| f A(P ) cm 

Here we have evaluated the addition to the tt — n scattering phase shift coming from the 
resonant scattering into the K meson state at a center-of-mass energy which corresponds to 
the laboratory energy Pq + APq determined by Eqs. 



dJP 2 -P 2 

A(P ) cm = -54- AP = 7AP0 = ± 7 |M| (54) 

Orn 
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K 



TT TT 
FIG. 1: The contribution of Hw to ir — ir scattering involving resonant production of a K meson. 
Because of the singular i^-meson propagator, the amplitude corresponding to this graph will be 
first order in Hw when evaluated at the center-of-mass energies m# ± j\M\ as required for our 
application. 

The infinite volume decay amplitude A appearing in Eq. |53] is normalized following the 
conventions of Lellouch and Liischer so that the corresponding decay width is given by: 

h 



r 



K- 



A 



(55) 



Finally, these results are used to evaluate the terms in the quantization condition, Eq. 0B1 
for the original finite volume P ^ 0, it — n system. The terms in this equation which are 
first order in the matrix elements of Hw are: 



-Ak r 



f d5 (k) 
l Ok 



+ 



h 141 



32irm 2 K A(P c 



cm/0 



Ak cm i — X 

I OK J k=k c 



(56) 



where relation between Ak cm and APq is determined by k cm = ^JP 2 /A — m 2 which implies: 



Ak r . 



PqAPq 
4k 



(57) 



Combining Eqs. 1521 and 1571 then gives the desired connection between A and 
M = (7i + 7i-\H w \K°): 



P \M\fd5 (k)\ 



± 



k \A\ 2 
"'cm I I 



P^M\ rd^y\ 

327rm^7|M| ~~ Ak cm I dk J k=k c 



(58) 



Ak cm I dk J k=k c 

If this equation is solved for \A\ 2 we obtain the desired generalization of the original Lellouch- 
Liischer condition to the case of P ^ 0: 



dq 



(59) 



This formula differs from that of Ref. |2j by the presence of the factor of 7 2 and the more 
complex, 7-dependent definition of the function (f>(q) given in Eqs. ETIandl 
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V. CONCLUSION 



In the preceding sections we have examined the case of two interacting particles confined 
in a finite spatial volume and carrying non-zero total momentum. We have determined a 
relation between the quantized energies of these finite-volume states and the two-particle 
scattering phase shifts. This result, first obtained by Rummukainen and Gottlieb in Ref. 
is here obtained from the Bethe-Salpeter equation of relativistic quantum field theory using 
an extension of the methods that Luscher applied to the case of zero total momentum in 
Ref, QQ. We then exp.oit these finite vo.n m e r e 8 „.t 8 to a „a. yZ e two-part,* deeays. The 
result, an extension of earlier work of Lellouch and Luscher, Ref. Q], to the case of non- 
zero total momentum, provides an explicit formula that relates finite volume decay matrix 
elements computed using lattice gauge theory techniques and the infinite-volume quantities 
that enter physical decay rates. 

The ability to work with states with non-zero total momentum when computing such 
decay matrix elements offers two potentially important benefits to the study of K — > im 
decay. First, by a proper choice of the total momentum, the corresponding two-pion state 
with lowest energy can be arranged to have an energy equal to that of the K meson, permit- 
ting the direct calculation of a physical, on-shell, decay matrix element. Second, by working 
with a K meson with non-zero momentum, we insure that the unphysical vacuum decay 
amplitude, normally dominant in such a Euclidean matrix element calculation, will vanish 
because of momentum conservation. We are now exploring this approach numerically. 

APPENDIX A: SINGULAR PART OF TWO-PROPAGATOR PRODUCT 

The singular contribution to the integral in Eq. will come as the energy Po is adjusted 
to cause singularities present in both of the single-particle propagator factors, A((P/2±/c) 2 ), 
to pinch the ko contour. Thus, we can obtain the singular part of the integral over ko in 
Eq. EH] by replacing both propagators by their singular part: 




(Al) 
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Equation |2U is based on the following expression for the singular part of the product of two 
free scalar propagators: 

* i -in6(P-k) 

(P/2 + k) 2 -m 2 + ie (P/2 - k) 2 —m 2 + ie P 2 /4 + k 2 - m 2 + ie 1 ' 

where these two quantities are equivalent in the sense that they have the same pole in the 
total energy P at the point Pq = u + + oj- 

Following Luscher, this formula is interpreted as relating two distributions in the variable 
k Q over the space of test functions f(ko) analytic in k within a band around the imaginary fco 
axis. Their equivalence can be demonstrated by multiplying by such a test function and then 
integrating fc along the imaginary axis. The left-hand-side of Eq. lA2l is evaluated by moving 
the k contour past one of the two pinching poles, k = P /2 — u_ or k = —Po/2 + uj + and 
keeping the contribution of that pole given by Cauchy's theorem. Following this procedure 
the left-hand-side of Eq. IA2I becomes: 

2u + uj_{P — oj + — uj_) 

where we have simplified this expression by replacing Pq in the residue of the Pq = u + + co>_ 
pole by its value at that pole. Next we evaluate the right-hand-side of Eq. IA2l bv using the 
delta-function to evaluate the integral over k . We obtain: 

-iTr — p ) (A4) 



P {p 2 /4+(^) 2 -P 2 /4-k 2 



m 2 + ie 



where for clarity the Pq arguments have not been simplified. The final step requires comput- 
ing the location and residue of the pole in Pq found in the expression in Eq. IA4I The result 
agrees precisely with that shown in Eq. IA3I provided one is careful to include the effects of 
the (P ■ k/Po) 2 term when computing the residue and then replaces Pq in the residue with 
u + + u_, it value at the pole. We forego a discussion of the specific region of analyticity for 
the test functions above and its consistency with the analytic properties of the amplitudes 
in the Bethe-Salpeter equation, Eq. [THl since these kinematics should be a direct Lorentz 
transform of those discussed by Luscher in Ref. "\. 
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